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WELL-POSEDNESS FOR COMPRESSIBLE MHD SYSTEM WITH 
HIGHLY OSCILLATING INITIAL DATA 

JUNXIONG JIA, JIGEN PENG, AND JINGHUAI GAO 


Abstract. In this paper, we transform compressible MHD system written 
in Euler coordinate to Lagrange coordinate in critical Besov space. Then we 
construct unique local solutions for compressible MHD system. Our results 
improve the range of Lebesgue exponent in Besov space from [2, N) to [2, 2N) 
with N stands for dimension. In addition, we give a lower bound for the maxi¬ 
mal existence time which is important for our construction of global solutions. 
Based on the local solution, we obtain a unique global solution with high os¬ 
cillating initial velocity and density by using effective viscous flux and Hoff’s 
energy methods to explore the structure of compressible MHD system. 


1. Introduction 

Magnetohydrodynamics (MHD) is concerned with the interaction between fluid 
flow and magnetic field. The governing equations of isentropic compressible MHD 
system has the form 

dtp + div(pM) = 0, 

dt{pu) + div{pu 0 u) + VP(p) 

^ H • VH — iV(|Hp) +/lAu + (A +/r)Vdivw, 

dtB + {divu)B + u ■ VB — B ■ Vit = vAB, 
divB = 0, 

{p,u,B)\t=o = {po,uq,Bq), div Ho = 0, 

where p = p{t,x) denotes the density, x S t > 0, u G is the velocity 
of the flow, and B G stands for the magnetic field. The constants p and 
A are the shear and bulk viscosity coefficients of the flow, respectively, satisfying 
p > 0 and A + 2/x > 0, the constant u > 0 is the magnetic diffusivity acting as 
a magnetic diffusion coefficient of the magnetic field. P{p) is the scalar pressure 
function satisfying P'{p) > 0 where p is the equilibrium state of the density. 

MHD system has been studied by many authors for its physical importance and 
mathematical challenges. X. Hu and D. Wang constructed the global weak solutions 
with large initial data in a series of papers [ziiHiin]. In addition, they also studied 
the important low march number limit problem in their papers. In 2012, A. Suen 
and D. Hoff [T3] extends small energy weak solution theory about Navier-Stokes 
system to MHD case. Then global weak solutions with initial data contain vacuum 
also be constructed by many authors in [laiiH]. 
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Besides the studies about weak solutions mentioned above, there are also a lot 
of investigations about classical and strong solutions. In 1984, S. Kawashima stud¬ 
ied a very general symmetric hyperbolic-parabolic system which incorporate MHD 
system as a special one in his PhD thesis [28]. He constructed the global solution for 
MHD system when the initial data is small in After S. Kawashima’s inves¬ 

tigation, recently, many authors [201 EB Ea [23] studied the global well-posedness 
and optimal time decay rate for the MHD system with the magnetic equilibrium 
state to be zero. Their studies get the optimal time decay rate for the derivatives of 
the solution which generalize the optimal time decay rate results for Navier-Stokes 
equations to the more complex MHD equations. 

R. Danchin did a seminal work in [24] in 2000. By scaling analysis R. Danchin 
gave the definitions about critical space for compressible Navier-Stokes equations, 
then under this new framework he constructed a global solution when the initial 
data to be small in the critical Besov space. After this important work, there are 
a lot of further studies mini El]. Especially, C. Wang, W. Wang and Z. Zhang |6] 
construct a global solution for isentropic Navier-Stokes system which allow highly 
oscillating initial density and velocity in 2014. In [261EI]! D. Bian, B. Yuan and 
B. Guo, firstly, constructed the local solution for compressible MHD equations in 
critical Besov space framework. In m, C. Hao analyzed the linearized hyperbolic- 
parabolic system related to MHD equations carefully, then he constructed the global 
solution in critical Besov space framework. Recently, there are some important 
works for incompressible MHD equations [101 [min]- In these excellent papers, 
F. Lin, L. Xu and P. Zhang constructed global solutions for incompressible MHD 
equations without magnetic diffusion by cutely using the dissipative structure of 
the system and the anisotropic Besov space technique. Then, X. Hu in paper |13] 
constructed the global solutions for compressible MHD equations without magnetic 
diffusion. Due to the structure is very different for incompressible and compressible 
MHD equations, X. Hu’s work is very different to F. Lin, L. Xu and P. Zhang’s 
work and X. Hu used lots of structure information about the compressible MHD 
equations. 

In this manuscript, we plan to construct global solution for MHD system with 
highly oscillating initial velocity and density. For isentropic Navier-Stokes system, 
global solutions with highly oscillating initial velocity can be obtained by a careful 
analysis about Green’s matrix for the linear system mm. In paper [33] , the authors 
studied the global well-posedness for compressible viscoelastic fluids in based 
critical Besov space by decomposing the whole system into three small systems with 
each one similar to Navier-Stokes case. However, as mentioned by A. Suen and D. 
Hoff in their paper [14], the presence of magnetic field effects results in a different 
and more intricate coupling of rates of partial regularization in the initial layer, 
so we need more careful analysis, we firstly transform compressible MHD system 
written in Euler coordinate into Lagrange coordinate in critical Besov space which 
is not appeared before. Then based on Lagrange coordinate form, we construct 
the local solution. Gomparing to the known local results [26l El] , we improve the 
uniqueness range for p from [2, N) to [2, 2N) where N stands for dimension and 
we also give a lower bound for the maximal existence time which is important for 
our global result. Based on the local well-posedness, we obtain a global solution by 
using effective viscous flux and Hoff’s energy estimates to explore the structure of 
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the system. For more explicit statements, please see Theorem 12.61 and Theorem l2.8l 
in next section. 

The other parts of this paper are organized as follows. In section 2, we give a 
brief introduction about Besov space and the two main results of this paper. In 
Section 3, we prove the local well-posedness in critical Besov space framework by 
using Lagrange coordinate methods and prove that the local solution can propagate 
the smoothness of the initial data. In Section 4, we use Hoff’s energy methods to get 
a uniform estimates about the solution. In Section 5, we prove a blow up criterion 
and then combine the results from Section 3 to complete the proof of Theorem 
12.81 At last, for the reader’s convenience, we collect a lot of useful Lemmas in the 
Appendix. 


2. Main Results and Some Preliminaries 

In this section, we will introduce some notations and give the three main results 
of this paper. Firstly, let me give some basic knowledge about Besov space, which 
can be found in [T] . The homogeneous Littlewood-Paley decomposition relies upon 
a dyadic partition of unity. We can use for instance any (j) G supported 

in C := {^ G 3/4 < |^| < 8/3} such that 

^</>(2-«C) = l if 

q& 

Denote h = T~^4>, we then define the dyadic blocks as follows 

AqU := (j){2~'^D)u = f h{2'^y)u{x — y) dy, and SqU = AfcU. 

The formal decomposition 

U = AqU 
gSZ 

is called homogeneous Littlewood-Paley decomposition. The above dyadic decom¬ 
position has nice properties of quasi-orthogonality: with our choice of (/, we have 

Aj^AqU = 0 if |A: — gl > 2, 

Ak{Sq^iAqu) = Q if \k-q\>5. 

Let us now introduce the homogeneous Besov space. 

Definition 2.1. We denote by 5/ the space of tempered distributions u such that 

lim SjU = 0 in S'. 

j-y-oo 

Definition 2.2. Let s be a real number and (p, r) be in [1, oo]^. The homogeneous 
Besov space R®,, consists of distributions u in 5/ such that 

<+ 00 . 

From now on, the notation R® will stand for R® ^ and the notation R® will stand 
forRIi. 
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The study of non stationary PDE’s usually requires spaces of type L^{X) '■= 
L^{0,T;X) for appropriate Banach spaces X. In our case, we expect X to be a 
Besov spaces, so that it is natural to localize the equations through Littlewood- 
Paley decomposition. We then get estimates for each dyadic block and perform 
integration in time. However, in doing so, we obtain bounds in spaces which are 
not of type This approach was initiated in [5] and naturally leads to 

the following definitions: 

Definition 2.3. Let {r,p) € [l,+oo]^, T G (0,+oo] and s G R. We set 

ll“ll4(A) 



and 


LJriB;) -.= {ueLJr{B;),\\u\\ 




< +00 I 


Owing to Minkowski inequality, we have Llj.{Bp) ^ Ll^{Bp). That embedding is 
strict in general if r > 1. We will denote by CriSp) the set of function u belonging 
to L^{B;)nc{[o,nB;). 


We will often use the following interpolation inequality: 








with 

^ Q ^ _ Q 

- =-1-and s = 9si + (1 — d)s 2 , 

r ri r 2 

and the following embeddings 

-A L^(Co) and -A C'([0, T] x R^). 

Another important space is the hybrid Besov space, we give the definitions and 
collect some properties. 


Definition 2.4. let s,t G R. We set 

and 

(R^) := {u G < Too} , 

where Rq is a fixed constant. 

Lemma 2.5. 1) We have B^’^ = B‘^ ■ 

2) If s<t then = B^p O y. Otherwise, Bf,’l = B^p + B^p. 

3) The space Bp’p coincide with the usual inhomogeneous Besov space. 

4) If Si < S 2 and ti > <2 then BffJ^ ^ 

5) Interpolation: For si, 82 , 01,02 G R and 9 G [0, f], we have 

ll/ll D»U + (l-f>)02.«<^l+(l-8)<^2 < ll/llsn.'^l ||/||^S2.<'2 • 

^2,p ^2,p 
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From now on, the notation S®’* will stand for B®’* and the notation i?®’* will 
stand for 

Throughout the paper, we shall use some paradifferential calculus. It is a nice 
way to define a generalized product between distributions which is continuous in 
functional spaces where the usual product does not make sense. The paraproduct 
between u and v is defined by 

TuV ■— ^ Sq-luAqV. 
qSZ 

We thus have the following formal decomposition: 

uv = TuV + TyU + R{u, v), 

with 

R{u,v) := '^AquAqV, Aq := Aq_i + A, + A^+i. 

q& 

We will sometimes use the notation T^v := T„u + R{u, v). 

For more information about Besov space and hybrid Besov space, we give refer¬ 
ence mElSE]. Throughout this paper, C stands for a “harmless” constant, and 
we sometimes use the notation A < B as an equivalent of A < CB. The notation 
A» B means that A < B and B < A. 

With these preparations, we can state our two main results. 

Theorem 2.6. Let p > 0 and cq > 0. Assume that the initial data {po,uo, Bq) 
satisfies 

Po- pG co< po< Cp \ 

Mo, Bo S B^ -^ 

Then there exists a positive time T > 0 such that if p G [2,27V), the system has a 
unique solution {p — p, u, B) satisfies 

p-pG C([0, T]-,B^/,^), ico < p < 2co-\ 

u,Bg C{[0,T];B^//-^)nL\0,T-B^//+^). 

In addition, we can take the maximal existence time T as follows 


(1 + ll“0|l5N/p)4’ 

P>1 

where c is a small enough positive constant. 

Remark 2.7. The above results are necessary for our following global well-posedness 
results, however. Theorem 12.61 has its own interest. As C. Noboru and R. Danchin 
[29| improved the uniqueness condition for Navier-Stokes system, we improve the 
uniqueness condition from p G [2, A) to p G [2, 2A) by using Lagrange 

coordinate technique. In our case, we need to get the magnetic field equation in 
Lagrange coordinate and estimate more complex coupling terms. In addition, we 
get a lower bound for the maximal existence time which is necessary for our global 
existence theorem. 
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Theorem 2.8. Suppose the equilibrium state p > 0, B = 0 and dimension N = 3. 
Let Co to be a small constant, 

P+ = sup |P^^^(/ 9 )| and P_ = inf \P'{p)\. 

co/ 4 <p< 4 c-i '=o/ 4 <p< 4 c-^ 


Assume that the initial data (pq,uo) satisfies 

po-pG Co < po < 

uo,BoGH^-^nBl^f-\ 

for some p G (3, 6 ) and s > 3. There exist five constants ci, C 2 , C 3 such that 

C 2 ^3 

IIpo - pIIl 2 < Cl, ||uo||g3/p-i < , I,—T]- \\Bo\\g3/v-i < Tj——n—fi- 

^P,i (1 + ||ao||^3/p)5 ^p,i (1+||ao||g3/p)2 




\Bo\\jj-s < 


C2 


(1 + llaoll^3 /p)5(1 + IIpo - pWI. + \\uo\\%^ + \\Bo\\%^) 

pA 

_C 3 _ 

(1 + llaoll p 3 /p)^(l + IIpo — p|||i2 + ||mo|| + llaoll 1 ^ 2 ) 


for some S G “ |i f); then there exist a unique global solution (p, u, B) satisfying 

p>^, p-pG C'([ 0 ,oo);i 7 ®), 
u,Bg Ci[0,oo); H^-^) n L‘^{0,T-, H^), 


for any T > 0. 


The proof of the above theorem used the advantage of harmonic analysis and 
the structure variable called effective viscous flux. The magnetic field and velocity 
field coupled with each other in the compressible MHD system so we need to es¬ 
timate many coupled terms. Based on the Theorem 12.61 we used energy methods 
in |14] which can incorporate the structure conditions of our system. However, our 
regularity conditions are better than the conditions in |14] , so we need not add ad¬ 
ditional conditions for viscosity coefficient to extend the local solution to a global 
one. 


Remark 2.9. Isentropic Navier-Stokes system is well-posedness when the initial 
velocity field have large oscillation which is proved in HIS]. Until now, there are no 
such results for compressible MHD systems for the coupling between velocity and 
magnetic field make MHD system are more complex than Navier-Stokes equations. 
The above theorem allows the initial velocity field and magnetic field have large 
oscillation and allows the initial density to have large oscillation on the set of small 
measure. Comparing to previous results for MHD system, the theorem above highly 
reduce the requirements about the small condition on the initial data. 


3. Local Existence and Propagation of Regularity 

In this section, we will prove Theorem 12.61 and prove that the obtained solution 
can propagate the smoothness of the initial data. Without loss of generality, we 
can assume p = 1 from this section to the last section. 
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Let Xu be the flow associated to the vector held u, that is the solution to 

( 3 . 1 ) Xu{t,y) =y-\- [ u{T,Xu{r,y))dT. 

Jo 

Denoting 

Kt,y) = u{t,y) = u{t,Xu{t,y)), 13{t,y) = B{t,Xu{t,y)). 

Then we reformulate the MHD system in the Lagrangian coordinate to obtain 

dt{JuP) = 0, 

JuPodtU — div {ad]{DXu) {2fiDA^u + Adiv^^uld + P{p)ld)) 

= div (^adi{DXu)BB'^^ - ^div (^ad]{DXu)\B\‘^'j , 

JudtB - udiv (adj{DXu)AlXB^ 

= div (^adi{DXu)Bu^^ — div (^adi{DXu)uB'^ , 

div{AuB) = 0, 

where Au = {DyXu)~^. During the transformation between Euler coordinates and 
Lagrangian coordinate, we used (ED). 

Next, we introduce ul the “free solution” to 

J dtUL - p-Aul - (A + p)VdivuL = 0, 

|^ML|t=0 = Uo, 

Tree solution” to 

r dtBL - i^ABl = 0, 

< divHi = 0, 

[^L|t=o = Bq. 


Lpg{u) = dtu -div {2pD{u) + Adivuld), 

Po 


Lu{B) =dtB-vAB. 

We can rewrite the second and third equations in system (IX^ as follows 

Lpo{u) = /i(u,m) + Po Miv(/2(M,M) + h{u,u) + h{u)) 

+Po Miv(J5(M, B) + Iq{u, H)), 
Lu{B) = Iriu, B) + div(J8{u, B) + Iq{u, B) 

+Jio(m, B, u) + Iii{u, B, u)), 



(3.3) 

and Bl the 

(3.4) 


Denoting 

(3.5) 

and 

(3.6) 
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where 


Ii{v,w) = (1 - Jv)dtw, 

l 2 {v,w) = (adj(-DX^) — Id){fiDw ■ Ay + fiAy\7w + XAy : Vw), 
I^iy, w) = fiDw ■ {Ay — Id) + ii{Ai^ — Id)Vw + X{dL^ — Id) : Vw, 
h(v) = a.d]{DXy)P(Jy^po) Id, h{v,w) = a.d]{DXy)ww'^, 

h{v,w) = ]^a.<li{DXy)\w\^, h{v,w) = (1 - Jv)dtw, 

Is{v,w) = v{a,d]{DXy) — Id)AyVw, Ig{v,w) = — Id)Vw, 

Iio{v,w, k) = a.d]{yXy)wk'^, Iii{v,w,k) = a,d){DXy)wk. 


In order to prove the local well-posedness, we define a map $ : (w, b) i-A {u, B) 
where (u, B) satisfies (13.711 with the u, B on the right hand side changed to v, b 
separately. Denoting u = u — ui, and B — B — Bj^, then we know that u and B 
satisfies 


(3.8) 


Lpo{u) = {Li - Lpg)uL + Ii{v,v) 

+Pd^div{l2{v,v) + l3{v,v) + hiv) + l5{v,b) + lQ{v,b)), 
< L,j{B) = Ir{v,b) +div{Is{v,b) + Ig{v,b) + Iio{v,b,v) + Iii{v,b,v)), 
div(.B) = div((/d — Ay)b), 
iu,B)\t=o = i 0 , 0 ). 


In order to give a clear statement, we need to define 

Ep{T) = [ve L{[0,T]-,B^/,P-^) : dtV,V\ G L\0,T-B^^{^-^)] , 
Bep(t){ul,Bl-,R) = {v,b: ||?; - UL\\Ey(T) + \\b- Bl\\ep{t) < • 


Suppose that T,R< 1, and T, R small enough such that 

rT 


(3.9) 


f \\Vv\\^N/p dt < c < 1. 

Jo p-i 


Step 1 : Stability of Bg (x)(ul,Bl;R) for small enough R and T. Us¬ 
ing Lemma mu we could get 


||u||ep(t) <Ce'^<’0''^'^(^\\Ii{v,v)\\^i^^^N/^p-i^ + ||(Li - 

+ (1 + II«o|IbW/p)( 11/2(1;,-f 11^3(1;,f)|li^(BW/pj 

Using Lemma [6.21 we could get 
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There are many terms can be estimated similar to Navier-Stokes system as follows 

WiLi- Lpo)ul\\, gN/p-i. < C'||oo||giV/p(l + ||ao||^iV/p)||'UL||^1 ,gN/p+is, 

T"'• p,l ' p,l Pil T ' p,l ■' 

||/i(ii,w)||^^(^N/P) < C{1 + l|ao|lB«/p)lklli^(B«/r)lklli^(B^/P), (i = 2,3), 

< CT{1 + llooll b))^/p)(1 + 

Here, we only give the estimates for the terms different to Navier-Stokes systems. 
l|4(v,w)|li^(B^/P) < C (||adj(i:>X^) - -h l) llw|||^(gN/p^ 

< (7(1 -I- lkllL:^(B"(P+7)ll^llL|(B"(P)’ 

||76(r’,w)||^^(3N/p^ < <^(1 + ll^^lll,^(B"/P+i))lkll||(B^/P)> 
\\I7{v,w)\\li^^^n/^p-i^ < (7||1 - -^p|||,~(B^/P)l|9tw|l^^(^N/p-i) 

ll'^8('i’, 11^)11 — C^lladj(ZlTCp) /(i||ll^j, ll^^ll 

+ (7||adj(I?Xp) — 

||/9(z;,w;)||^^(^«/P) < C'Pp --fcl|||,^(B^^/p)lkllB^(B"/*>+7 


<C^IHI 




llIio(v, w,k)H^,^^j^N/p^ < (7(1 -k ||adj(V7C„) -/d||^^^^^/pp||w||^^^^^/p^||fc||^^^^^/p^ 
< (7(1 -I- lkllB;^(BW/P))PllB2^(BN/p^||'U;||£^|.^N/pj, 

Combining the estimates about Ii to In, we could finally get 

||m||bp(t) <(7e'^'’0-’”'^(l -I- ||ao||^N/p)^|T -|- ||ao||^N/pH ulII^^^^n/p+i^ 

and 


\\B\\e^{t) <C'(^(||ml||^^j.^n/p+i^ 

where ii = v — ul, b = b — Be- 
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We could reduce the above two inequalities further to get 
||■w||_Ep(T) + ||ao||^N/p)^|T + ||ao|| 

(3.10) +iR+ 


and 


(3.11) 


\\B\\Ep(T) <C{R+ l|wL||^|(gN/P) + ll“i|lL;^(B"/P))(l^ + 


Here, if we assume 

Cpg^mT < log2, T < R^, 

W^oW b^/^p\\ul\\ ii^0N/^p^ < R^, 

W^tRLWll^^gN/^P-l) + II^l||£^(b"/p+1) + \\Bl\\i2^^bN^ < R, 

(1 + WaoW^N/^pfR < 77 < 

then we know that 


I|m||bp(t) < R, I|i?IIbp(t) < R- 

Hence, $ is a self map on the ball Be^{t){ul,Bl',R)- 

Step 2 : Contraction estimates. Taking (t;^,6^) G Bep(t)(ul,Bl;R) and (v^,b^) G 
Bep(t){ul, Be': R), we define {u^,B^) = ^{v^,b^) and (u^, B^) = ^{v^,b^). Denote 
Su = u'^ — , SB = B'^ — B^, Sv = v'^ — and Sb = b‘^ — b^. 

Through simple calculations, we have 

LpoiSu) =Ii{v^,Sv) + {Jyi - Jy2)dtv‘^ 

+ Po ^{(d2(^’^^’^) - l2{v'^,v'^)) + (/3(^’^^’^) - /3(u\u^)) 

+ ihiv^) - hiv^)) + {h{v\ b^) - h{v\b^)) 

+ (/6(^'^fe^) - /6(n\5i))|, 

and 

L^{SB) =(17(7;^, 6^) _+div|(/8(n^,6^) -/ 8 (n\ 5 ^)) 

+ ihiv^^b"^) - Ig{v^,b^)) + {hoiv^, b'^, n^) - Iio{v^,b^,v^)) 


+ {Iii{v‘^,b'^,v‘^) - Iii{v^,b^,v^ 
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Under the condition Cp^^mT < log2, using Lemma [Q] and Lemma [6?2l we get 

+ ||(J„1 - + \\l 2 {v^,v^) - h{v'^ 

+ - l3iv^,v^)\\i^^gN/p^ + \\hiv^) - hiv^)\\i^^gN/^p^ 

+ ||/5(p^6^)-/5(P^6^)|li^(^./P) 

+ Whiv^, b"^) - 


and 


\\SB\\E^iT) <C[\\lAv\b^) - l7iv\b^)\\ 

+ \\hiv^,b^) - h{v\b^)\\ii^f^gN/^p^ 

+ \\I9{V^ ,b^) - l9{v\b^)\\li^^BN/^P) 

+ \\hi{v‘^,b^,v^) - 

There are many terms appeared in the Navier-Stokes system, we only list the esti¬ 
mates 

\\Ii{v^,Sv)\\^i^^^N^p-i^ <C\\v^\\^^^^N^p+i^\\dtSv\\^i^^^N/p-iy 

||(J^i Jy2)dtV ^C||<^'Il||^^^^N'/p+i^ ||(3tP I 

\\l2iv^,v^) - l2iv\v^)\\^,^^^N/^p^ <C||(p\i;2)||^^(^«/p+i^||(5i;||^^(^«/p+i^, 

\\l 3 {v^,v^) - l 3 iv\v^)\\ <C||(p\ i;^)||||(5i;||^^(^^/p+i^, 

Whiv^) - hiv^)\\ii^(^BN/^p^ <C{1 + ||ao||^N/p)T||(5i;||^^^^N/p+ij. 

Now, we analyze the different terms carefully. Since 

hiv\b^) - hiv\b^) =(adjpX„2) - adi{DX,i))b^bY 

+ adj(i^X,i)(&2-6i)(&2)T 

+ adi{DX^i)b^{b^ - b^f, 

we have 

\\hiv'^,b'^) - h{v^,b^)\\ ^C'||(5u|| l:^(b^/p+1)II^^II|2^(s«‘/p) 

+C'(l + lk^ll£,^(B^/p+i))ll<5^lli2^(B«/P)ll(^\&^)lli2^(B"/P)- 
Iq can be estimated same as I 3 . Since 

/ 7 ( w ^, 6 ^) - l7{v^,b^) = {Jyi - Jy2)dtb'^ + (1 - Jv^){dtb'^ - dtb^), 
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we know that 

.,2 l 2 \ 


\\I 7 iv ,b ) - hiv ,h <C\\Jy2 - Jyi\\^^^^N/p^\\dtb 


+ CII J„i - l||£oo(B^^/P)||Gtoe||^^(^«/p-i) 


After some calculations, we have 

h{v^y) - h{v^y) =J^(adj(DX^2) - adj(£>X„i))A^2V&^ 

+ v{a.A]{DX^i) - Id){A'^2 - A^i)V6^ 

+ i2(adj(i:»X„i) - Id)A^^V{b^ - &i), 

so we can get 

Due to 

/9(u2,62) _ I^{v\b^) = iy{A^. - A^i)V62 + v{Al, - Id)VSb, 

we have 

||/9(i; 2,62)-/9(i;i,6^)11 <^II^^'II4 (b™)II^"II4(b"/-) 

Since 

/io(w^, 62,^2) - Iio{v^,b^,v^) =(adj(VX„ 2 ) - adj(VX^i))62(i;2)'^ 

+ adj(VX„i)(62-6i)(i;2)^ 

+ adj(VA:„i)6^(i;2 

we easily have 

\\hoiv‘^,b‘^,v‘^) - ^loiv^,b^,v^)\\ 

<C'||6u||^^(^n/p+1)||62||~^(^n/p)||v^||~^(^n/p) 

+ C'(l + ||v lli;^(B«/p+i))ll<^^lli2^(B"/p)lk lli2^(B"/P) 

+ (7(1 + l|v^|l£;^(S«/P+l))ll^^ll|,|(B^/P)l|dv||^^(^W/p^. 

The term ||/ii(w2^ 62, i;2) — 6^, v^llUi /rN/p, can be estimated similar to the 

^TV^P,1 

above terms, so we not write the details. Combining all the above estimates, we 
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finally arrive at 

||<5w||_Ep(T) <^(1 + ||ao ||{T 
+ \\dtv^ 




^AT/p+l. 




11'^ 

and 

WSBU^iT) <c{(||9, 




}' 


+11^' iIl^(b 






,N/p + l 


) 


+ II^^IIl|,(b"{'’)II'^^IIl|,(b"(p) + II^^IIl:^(b"('’ 

”*" 11'^ IIb^(b"{’’+^)II'^^IIl:^(b"(’’+^) “*" 11^ llB|(B^{^)l^'^^llB|(Bp;yj J 

Through some simple calculations, using conditions (13.121) with may be lager con¬ 
stant C, we have 

II'^“IIbp(t) < C{1 + ||ao||^iv/p)^i?(||(5i;||£;p(T) -I- ||^6 ||ep(t))i 

P ,1 

and 


cjAf/p+i^ ^^iv/p-ij 


,iv/p-l«^ 

,iV/p. 


\\^B\\e^(t) < CR{\\6v\\ e^(t) + II<56||bp(t))- 

Hence, we get 

\\{Su,6B)\\e^(t) < C{1 + ||ao||^w/p)^i?||((5u,(56)||Ep(T). 

P ,1 

Due to conditions (13.121) . we finally get 

(3-13) IK^U, (5H)||£;p(T) < 2ll('^^’'^^)llsp(T)- 

Now, by contraction mapping theorem, we know that $ admits a unique fixed point 
in BepBe; R)- 

Step 3 : Regularity of the density. Denoting a = p — 1, and we already 
know that p = J~^pQ. So we have 

n = i^u ^ ~ l)no + no + {Ju ^ “ !)■ 

Due to u G L^(0,T; we know that Ju^{t) — 1 G (^([O, T]; So we 

know that a belongs to (^([O, T]; If 

R< An^n\\ - 

4(1 -I- ||ao|| E^/p) 

p,i 

we have 


p{t) > (1 + no) — (1 + l|no||L“)(||Mz,||£^^^w/p^ -I- K) 

> (1 -I- Oo) — 2(1 -I- llooll AN/p)i? 

"p,i 

> 0 , 

where cq defined as in Theorem 12.81 and t G [0,T]. 

Step 4 : Uniqueness and continuity of the flow map. Consider two cou¬ 
ples of initial data (pj, Uq, Hq) and (pg, "noj ^o)- Define 5p = p^ — p^, 5u = — v} 
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and 6 B = , we can easily find the equations for Su and 5B. Using similar 

methods as Navier-Stokes system and estimate the terms containing B as in Step 
3, we can finish this part. For there are no new gradients, we omit the details. 

At this stage, by same reasons as in [2^, we can get the local well-posedness 
results for our MHD system. 

Step 5 : Lower bound of the existence time. Taking m to be a large con¬ 
stant fixed as before and (7 to be a large enough constant, let 

^ ^ Cpo, mT <log2, {l + \\ao\\gN/pfR<r], T<R^, 

(3.15) ^ 

ll«oll rN/pII^^oIUn/p-i < R^, ||uoIUn/p-i < R, II^oIUn/p-i < R, 

p,l p,l Pil PjI 

then the conditions (13.121) and (13.141) are satisfied. From the above conditions, we 
know that 


(3.16) 


T < 


(72(1 -I- ||ao|LN/p)'‘’ 


T < 


log 2 

COr, .m 


Next, we calculate Cp^^m carefully as follows 


Cpg^m ^ C 


( — 

Po 


dr 


AiV/p 


(3.17) 


<C [ ||5'^Vpo||^N/pdr 

Jo pa 


< CT 2 ^^\\ao\\l 


B. 


N/p 


Combining (13.161) and (I3.17|) . we know that we can take T as 


(1 -I- ||ao|| gK/p)'^ 

P>1 

where c is a small enough positive constant. At this point, we finished our proof of 
Theorem 12.61 

Next, let us go to the second part of this section to prove the solution in Theorem 
12.61 can propagate the smoothness of the initial data. 

Before our proof, we need to introduce some notations which can also be found 
in [61130]. Define a weight function {ujkit)}kez as follows 

a;fe(^) = ^2'=-^(l-e-"2''*)l/^ keZ, 

i>k 


where c is a positive constant. We easily know that for any k G h, 

uJk{t)<2, ujk{t) ^ uJk'it) iikr^k', 

(3-18) , 

Wfc(t) < 2^ ^ ujk'it) if fc > fc', 0 Jk{t) <ioJk'{t) iik<k'. 


Now, we can introduce the following weighted Besov space. 


Definition 3.1. Let s G R, 1 < p, r < -l-oo, 0 < T < -|-oo. The weighted Besov 
space Bp j.{oj) is defined by 

Rp,r{^) = ^ ■ ll/ll,.(1.,;) < +Oo| , 
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where 




:= \\{2’^‘^Wk{T)\\Akf\\Lp), 


Obviously, C Bp^^{oj) and 




< 211/11 


-B; 


We also need to define the Chemin-Lerner type weighted Besov spaces as follows. 

Definition 3.2. Let s £ R, 1 < p, g, r < +oo, 0 < T < +oo. The weighted 
functional space L^(Bp^^(uj)) is defined as the set of all the distributions / satisfying 

II/IIl^(B| ,,(0;)) •= II (2^'*Wfc(T)|| Afe/(t)||i<,(o_T;LP))fc||fr • 

After this short introduction about weighted Besov space, we deduce some prop¬ 
erties for the solution we obtained in the first part of this section. Obviously, we 
have 

(3.19) j<p(t,x)<2co\ 

Since 


C3 




and 


(1 + IkolIpN/p)^ 


^2 

Bp 

we know that 
(3.20) 


ll“llL5S(B-/'>-) - < (TTKI^ 

For the magnetic field, we have 

(3.21) IIRII AM/n-l- . + IIBII ;i , AAT/p+i, < i? -f llBnll aw/p < 


C2 + C3 


2c3 


+ ll^llLi,(Bf/'’+^) - R+\\Bo\\gN/p A ||ao||^iV/p)2' 

For the density, we know that 

II®IIl“(Bp"/0 - II^^IIb),(Bp"''oII“°IIbp"^'’ + llaolls^/P + II^w|Il^(b«/p)- 

Hence, by some simple calculations, we obtain 

(3.22) ll'^lliooc^^'/p., < 2||ao||^Ar/p -|- , |j 

) Bp (1-I-||ao||pN/p)^ 

So combining (13.191) . (I3.20L (13.211) and (13.221) . for dimension TV = 3, we know that 
the solution satisfies 

y < p{t,x) < 2co\ 

C 2 C 3 

II“IIl”(b^/ 0 - 2||ao|lB3/P + + 

(3.23) C2 C 3 

II^IIl”(b^/"- 0 + - (1 + ||ao|U 3 /p) 5 ’ 




mir. 


^ (l + ||ao||^3/p)2' 


2cp 
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At this stage, we estimate carefully. As in [ 6 ], we have 

^ 2^1||A,ao|Up<C2-5^||ao||/^., 

3>M 

and 

^ 2^-fc.,(T)||A,ao||LP < C22“r5||ao||ff2 +C2-5“||ao||ff2. 
j<M 

So if we take M large enough and 
(3.24) T = 


(l + ||ao||^3/p)4(l + ||ao|k2)i2’ 

we have 

(3.25) ||ao||^3/p(^) < c, 

where c, c' are small enough constants. Similar to (13.221) . we also have 

(3.26) 


C2 C 3 

^ ll«0lU3/r(^) + (l + ll^^ll )4' 


In the following proposition, we show that this solution allows us to propagate 
the regularity of the initial data in Sobolev space with low regularity. 

Proposition 3.3. Let p € (3,6) and 1 — | < d < |. Assume that {p,u,B) is a 
solution for system ()1.1|) . If {ao,UQ, Bq) £ x H~^ x then 

II“IIl“(S)-'’) - ^ (ll“o||ffi-« + lluollij-^ + ||i?o|lij-5) , 

(3.27) II('“>^)IIl“(b-2 ) 

< C ^I + ||ao||^3/p^ (11(^0, Bo)\\^-s + TP+||ao||//i-s). 

Proof. The proof can be divided into three steps. Set f) = 1 — S. Without loss of 
generality, we may assume that c, C 2 , C 3 , c', T (choose as in (13.251) and (13.241) 1 small 
enough such that 


C 


C2 + C 3 


(3.28) 


(1 + ||ao|U3/p)5 


C 


2 c 2 


+ 2 C c' + 


C- 


C2 + C3 


(I + ||ao||^3/p)4 
2 c 2 


(1 + Il«o||o3/p) (1 + ||ao|| A3/p) 


CTP+ < -, 


where C is the constant appearing in the following estimates. 

Step 1 : Estimate for the Transport Equation. This estimate is similar to 
the Navier-Stokes system, so we omit the details and only give the results as follows 


(3.29) 






For details, we give the reference [ 6 ] on page 192. 

Step 2 : Estimate for the Momentum Equation. Denote ft = - and A = 


p ’ 
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apply the operator Aj to the momentum equation of (HID, we can obtain 
dtAjU — div(/iVAjw) — V(AdivAjM) 

= AjG + div([Aj, /l] Vm) + V([Aj, AJdivrt) 


(3.30) 


where 


+ A, [-B S/B 


\ - Ia, 




-y\B[ 


G = —u ■ Vm — pVu + ^Vpdivu. 

p p^ 

Taking energy estimate for weighted Besov space, we have 


(3.31) 


Similar to the Navier-Stokes system, we have 

2^'%-(T) (||[A„/2]ViiL^(i.) + ||[A„/2]ViiL^(i.) 


+ c 


> 

■Ri 

<1 

^(L2) + II [A 

+ c 

-B-VB 

+ G 

-vd^n 


P 


p 


£2 






(3.32) 
and 

Then, we estimate the term which is not appeared in the Navier-Stokes system as 
follows 

^ B VB 


(3.34) 


Plugging estimates (I3.32p . (13.331) and (13.3411 into (I3.3ip . we obtain 
(3.35) + ^ll®llL“(B^/P(i.j))ll^llL:^(Bf_+^) + ^’^^+ll®llL“(Bf 2 ) 

+ <^(1 + ll«lli~(B|/r))||B||£2^(s3/p)||B||^2^(^3^). 

Taking energy estimates for equation (13.301) . we obtain 


(3.36) 


+ C 


2^^ (||[A„m]VwL^(,. 3) + ||[A„A]Vw|L^(i3) 


^2 


+ G 

-B-S/B 

+ G 

-vd^n 


P 


p 
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The first two terms can be estimated as in the Navier-Stokes system, we only give 
the results as follows 


(3.37) 


(3.38) 


+ ^-f+ll®llL“(B^ 2 )^' 


For the terms not appeared in the Navier-Stokes system, we can estimate as ()3.34p . 
Plugging p.37p . (13.381) and (I3.34p into (I3.36p . we have 

(3.39) + ^II'“IIl|,(b^^P)II“IIl|(b^,2 ) 

+ C'Cl + l|a|l|,op(B3/pp||S||^^^^3/p)||B||^^(^p^^. 

Step 3 : Estimate for the Magnetic Field Equation. Applying the operator 
Aj to the third equation of system m, we obtain 

dt^jB — i/AAjB = —Ajdiv{Bu^ — uB"^). 


Performing energy estimates, we could have 

(3.40) 


-®IIl“(b|-1) + II-®IIb^(b|+1) <CPo|Ib|-i 


Since 

(3.41) 


+ C||div(Bu -uB )||i^(B|-i). 




<C'||m||£2,(s3/p^||b||£2^(33^), 


and 

(3.42) 
we know that 

Combing estimates (13.351) and (I3.43p . we will get 

II^IIl^(B2^+1) + II«IIb|,(B2^_2) II^IIz),(B2^+1) 

— ^(Il^olls^-i + \\Bo\\gfi-i + TP+||a||£^^^3^^). 
Combining estimates (I3.39P and p.43l) . we will have 
ll^llL“(Bf-i) + \\^\\l^{B§-^) 


(3.43) 


(3.44) 


< Cdluoll^^-i -I- WBqW^p-i + 

+ C'||ao||^3/p||u||j^^(^3+i^. 


(3.45) 
























COMPRESSIBLE MHD SYSTEM 


19 


At this stage, combining (|3.44ll and (j3.45|) . we obtain 
(3.46) 




\B\ 




< (7(1 + llooll^s/p) {\\uo\\gP-i + WBoWgp-i + 

At last, combining (13.291) and (13.4611 . we finally get the desired results (13.2711 . □ 

In the last part of this section, we prove the solution propagate the regularity of 
the initial data in Sobolev space with high regularity. 

Proposition 3.4. Assume that {p,u,B) is a solution of system (11.111 on [0,T], 
which satisfies p > co, 

(3.47) aeZ??(B3/P), u,seZ??(B3/p)nL:^(B3/p+i), 

where a = p — 1. If (oq, mq, Bq) € x x for s > 3, then we have 

(3.48) aGL^iH^), u, B G (1 

Proof. Considering (I3.47L we can divide the time interval [0,T] into finitely many 
small intervals [Ti, Ti+i] with z = 0,1 ,..., IV such that 

ll(^) ^)\\L^(Ti,Ti+i;Bp^’’) — 

for some e small enough. Here w* = {w^} stands for 

< = ( 1 -e ) = . 

i>k 

Denote L^,{Bpp) := L'^(Ti,Ti+i] B^p). For the density, the estimate is similar to 
the Navier-Stokes equations, so we omit the details 

II®IIl|?.(B|_ 2) — C'(||a('Fi)||B|^^ + \\'‘^\\Lf,(B^+f-))- 

For the velocity field, as in the proof of Proposition 13.31 we have 

<C(||zz(r,)||^.-. + ||G||^. 


(3.51) 


C 


2^*a;, (t) ||[A„A]Vzz|U,^ + || [A,,/l]Vn|L^ 


£2 


+ c 

-B-WB 

+ C 

-V(|h2|) 


P 


P 




Similar to the Navier-Stokes system, we have 

(3.52) + + ll“llb5^(B|_2)ll^llL;^.(B|^^+^) 

+ (Ti+I — 7i)ll®llL“(B| 2))> 

and 


(3.53) 


2 ^-^a;j(t) (||[A,,A]Vn|L^^(i.) + ||[A„ 


£2 




(b: 


d,3/p + l 


,)■ 
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For the last two terms in 13.511 all can be estimated as follows 


1 


-B VB 


__ < C||B - VBII + Cll—B - 


(3.54) +c||B.VB||j._,j„,_.,Mj.|i.j ■ 

Plugging (I3.52L (13.531) and (13.541) into (I3.51|) . we will obtain 

II ^11 (4^+1) + ll^lli|,.(B|_2) 

< C[\\u{Ti)\\^B-i + ||a||£~(B|^^)(||w||(s3/P) + {Ti+i — Ti))) 


(3.55) 


+ C'||B||^^^|.^3/pj||B||£2^_(3.^^) + C\\ao\\^3 /p\\B\\ ^ 2 ^^^^yps^\\B\\ 

+ C||i?|||J(.3/p)||a||,.(^.^). 


For the equation of magnetic field, using similar methods as in Proposition 13.31 we 
obtain 


(3.56) 


Since 


ll'®llL“(B 2 ,’i^) —^II-^oIIbJ-i + C\\u ■ 

+ C||B.V«||£^(^p-4). 


11 ^■ '^■®IIl:^.(bj-1) -^(II^IIl2,.(b^''p)II^^IIl|,.(.b“-1) + II^^IIl|,.(b^''p-1)II^IIl^.(.bj 2 

<C'||m||£^^(^3/p)||B||^2^^(^. 2) + C'l|S||£^_^^3/p^ 11^11^2^^ 

and 

11^ ■ '^^IIl:^.(B^ 21) ^ ^(ll^llL|,.(B^/P)ll^lli|,.(B2“,2) ll“lli|,.(B^^*’)ll'®lli|,.(B|,2))’ 

we know that 

II^IIl“(B2»_-1) + II^IIl:^.(B 2"^') + ll-^llii.^iBi.s) 

(3.57) — ^ll'®(^i)llBj-i ^ll^llL|,.(.B|^'’)ll-®llh|,.(.B| 2 ) 

+ ^\\B\\ii^s^b1^i’)\\'^\\lI,(^.2)- 

Combining (I3.50F (13.551) and (I3.57F we have 

II«IIl“(B|, 2) + ll“lliF(S2n') + ll^lli|-,(4“.2) + ll^lliT..(4”^') l|■®lliT,(4^2) 


(3.58) 


< Ci\\am\sB+ ||w(r.)|lBrd + ll^(^OII( 


v^ 2,2 

-1). 


Estimate velocity as in (13.361) with regularity index to be s — 1 and combine the 
estimate about magnetic field B, we finally get 

ll'“llL|?.(B2,i^) ll'^lli?)(-B2,i^) -^(1 + II“IIl“(B^/P))(II“(^»)II.B|2 

+ ll«m)llBrd + ll^(^OIlBrd)- 


(3.59) 
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Combining (|3.58ll and (I3.59L we arrive at 

l“IIZ“(B|2) 11(^1 + ll(^’ 

< c(i + iiaiii=o(B3/.))(ii«m)iiB|_, + 

By induction, we know that 

Hence, we complete the proof. □ 


(3.60) 


4. Hoff’s Energy Method 


In two papers [mill], D. Hoff construct a global weak solution for Navier-Stokes 
equations with discontinuous initial data with small energy. In the paper m, A. 
Suen and D. Hoff generalize the results for Navier-Stokes equations to compressible 
MHD system. Here, we use the idea to our case. Comparing to mi, we remove the 
restriction on viscosity for we have stronger condition on the initial data. 

We set a{t) := min(I,t), define 


A,iT)= sup i\\Vumh + \\^Bit)\\l.) 

0<t<T 


10 dR3 


p\uf + \Btf dxdt, 


A 2 {T) = sup (cr(t) / p\u\'^ dx + a{t) [ \Bt\'^ dx] 

0<t<T V dR3 JR3 / 


10 ./R3 


cr(t)|Vu|^ dxdt ■ 


10 JR3 


cr(t)|Vi?t|^ dxdt, 


and 


E{T) = [ a{t) (IVHHhP + + lyuHsp) dx, 

H{T)= I I \VBfdxdt+ I I \Vufdxdt 
Jo dR3 Jo dR3 


-f 


cr(t)|VH|‘* dxdt 


cr(t)|Vu|‘* dxdt. 


In the following sections, we denote 

f = ft+u-Vf, p' = X + p, 


and 


Co = \\po-p\\h + \\uorH. + \\Borm. 

Throughout this section, we denote by C a constant depending only on A, p, cq, p, P+ 
and PZ^. 


Theorem 4.1. Let {p,u,B) be a solution of system U.l\) satisfying 

p-p€ ^([O, T]; iZ^), u,B€ ^([0, T]; iJ^) n L‘^{0, T; H^). 
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Then, there exists a constant eo depending only on v, X, fj,,co, p, P+, P_^ such that 
if the initial data {po,uo, Bq) satisfies 

Co < Poix) < Cp \ X e 


IIp-pII 


L2 


ffi + ll^ollffi < eO) 


then we have 


y < p(i, a;) < 2co \ {t,x)e[0,T]x 


A{T) := A,{T) + A^iT) < el 

Proof Considering the assumption, there exists a 0 < Tg < T such that the solution 
{p, u, B) satisfies 

^ < p{t,x) <2cl, (t,x) e [0,ro] X 

Ai{TlpA2{Tl<el 

Without loss of generality, we assume that Tq is a maximal time so that the above 
inequalities hold. In the following, we will give a refined estimates on [0, T] for the 
solution. Due to the proof is too long, we divide the proof into several lemmas. 

Lemma 4.2. (Lf energy estimate) 

[ \p-p\^ + pIu]"^+ \B\'^dx+ [ [ \Vu\^ + \yB\^dxdt<CCo. 

7r3 Jo Jr3 

The proof of this lemma is exactly the same as Lemma 2.2 in |19j . so we omit it. 
Lemma 4.3. (PI^ energy estimate) 

Alin) <CCo + C [ ° [ \yBf\Bf + \VBf\uf + \Vuf\Bf dxdt 
Jo Jr3 


rTo 


c 


|Vu|^ dxdt. 


Proof. Multiply the second equation in dm by ii and integrate over R^, we will 
obtain 


(4.1) 


/ p\u\^ dx = / (—ii • VP + +/r'Vdivui 

Jr^ Jr^ 


— V(i|Pp)ii + div(P B^) u) dx. 


By the continuity equation, we have 

(4.2) 9P + div(itP) = divu (P — P'p). 

There are some terms can be estimated as in the Navier-Stokes equations, so we 
omit the details and only give the results as follows 

(4.3) j —u-\7Pdx<dt I divM (P — P(p)) da; + C||Vm||| 2 , 

Jr3 Jr3 

(4.4) p ( Auiidx< — ^dt [ iVupdx + C f |Vupdx, 

7r3 2 Jr3 J^3 
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(4.5) /i' f ydWuiidx < —^dt f \divu\'^ dx + C ( \yu\^dx, 

Jr3 2 Jr3 J^3 


and 


(4.6) 


/R3 


divu [P - P{p))dx < ||divu||i 2 ||P(p) -P(p)||i 2 


<CQ\\divu\\L^. 

The following terms are not appeared in Navier-Stokes system 


(4.7) 


f — V(i|i3p) wda; + f div{B B^) iidx < C ( |Vi?||i3||u| da; 

Jk? 2 Jr3 

<C [ \VB\^\B\^dx + e [ lupda;. 

JR 3 Jr 3 


/R3 aR3 

Multiply the third equation in ( 11 . 11 ) with Bt and integrate over we have 

/ \Bt\'^ dx + ( div{Bu^ — uB^) Btdx =—^dt f |Vi?pdx. 

dR3 dR3 2 J^3 

Integrate the above equality from 0 to Tq, we will obtain 


(4.8) 


^ [ \VB\^dx+ [ “ [ \Bt\^dx 
^ dR3 Jo Jb .3 

< CCo +C [ ° [ \VB\^\u\^ + \Vu\^\B\^ dxdt. 
Jo Jr3 


Combining (14.IL (14.3L (14.41) . (14.5p . (j4.6p . (j4.7l) and (14.8L we finally get the desired 
result. □ 


Lemma 4.4. (H^ energy estimate) 

A 2 {To) <CCo + CAi{To) + C [ " [ a\Vu\^dxdt 

Jo Jr3 

+ C [ [ a{\u\^+ \B\^){\Bt\^+ \u\^+ \u\^{\yB\^+ \Vu\^))dxdt. 

Jo dR3 

Proof. Take material derivative to the second equation of dm to obtain 

put + pu-\7u + yPt + div(VP 0 zi) + V(i|P|^)t + div(M • V(i|P|^)) 

— div(P P^)t — div(M div(P B'^)) = pAut + pdiv{u Au) 
+ /i'Vdivuj + ^'div(u • Vdivit). 
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Multiply a{t)u on both sides of the above equality, we will have 


p|up dx 


— \<x' [ p\u\'^ dx = — I aii^djPt + div{djPu)) dx 

2 Jr3 Jr3 

= - f aum^\B\'^)^^^t+div{udji^\B\^)))dx 
+ f au^div{B^B)t + div{udiv{B^B))) dx 

Jr3 

+ pa / u (Am* + div(AM 0 46)) dx 

Jr3 

+ p'a / 46 (Vdiv46t + div(Vdiv46 0 u)) dx. 

Jr3 


Similar to the estimate, there are some terms same as Navier-Stokes equations 
and we only give the results as follows 

(4.10) — / aii^djPt + div{djPu)) dx < Ca\\Vu\\i^ 2 \\Vu\\i^ 2 , 

Jr3 


(4.11) pa f ii {Aut + div{Au (Si u)) dx < a f —^iVri^dx+ (tC j |V 46 |^dx, 

jR3 JrS 4 dR3 

pa / 46 (Vdiv46t + div(Vdiv46 0 u)) dx 

(4.12) f X u 

<a — — |div46p + —iVup + dx. 

Jr3 2 4 

For the terms not appeared in the Navier-Stokes equations, we have 
- / (^u^dj{^\B\%+diYiudj{hB\^)))dx 

(4.13) 

<e a\Vu\^dx + C alBl"^HBtl"^ + \uf\VBl"^) dx. 

JR 3 Jr3 

For the equation about magnetic held, multiplying aBt on both sides, we will have 


-cr / lUtP dx -I- 64 


/o dR 3 


a\VBt\'^ dxdt 


n l f 

aBt{div{Bu^ — uB^))tdxdt + - / / a'\Bt\^ dxdt. 

:3 2 Jq J^3 
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For the first term on the right hand side of the above equality, we have 

rTo 


(4.15) 


<e 


<e 


— f f aBt{div{Bu^ — uB"''))tdxdt 

Jo dR3 

/ f a\7Bt{Bt + B uj — Ut B'^ — u Bf ) dxdt 

Jo Jr3 

f f a\VBt\'^ dxdt + C ( f a{\Bt\'^\u\'^ + \B\'^\ut\'^) dxdt 

Jo dR3 Jo dR3 


rTo 


lO JR3 

+ c 


dxdt 
rTo r 


h dR3 


a{\Bt\M + \BV\W + \B\M\^uY)dxdt 


where e is a small enough positive number. Combining estimates (14.91) . (I4.10L 
(I4.11|) . (14.121) . (I4.13L (14.141) and (I4.15L we will obtain 

a f p\u\'^ + \Bt\'^ dx + f f a{\yu\'^ + \VBt\'^)dxdt 

Jr3 Jo dR3 


rTo 


<C 


(4.16) 


/O JR3 
/*1ATo 


a\yu\i^o\'\/u\i ^2 dxdt + C 


rTo 


10 JR3 


cr|VM|^ dxdt 


r*lATo 


C 

C 


/ / p\uf dxdt + C 

lo dR 3 Jo ^R 3 

rTo r 

|2 I Id|2\/'|d| 2 I I -12 I L,|2nv7D|2 , |V7„,|2 


Jo JR3 

We can easily obtain the following estimates 
rTo r 


Bt\ dxdt 

cr(|u|^ + |il|^)(|ilt|^ + \u\^ + |Mp(|Vil|^ + \Vu\^))dxdt. 


lo 


(4.17) 


a\Vu\i^ 2 \Vu\j ^2 dxdt 

rTo 


< 


C f f alVu\^ dxdt + C f ( a\Vu\’^ dxdt 
Jo Jr 3 Jo Jr 3 


< cc'o + CAi(ro), 


(4.18) 


pIATq 


/]R3 


p\u\^ dxdt - 


^IATq 


/1R3 


\Bt\‘^dxdt < CAi{To). 


Substitute (14.171) and (14.181) into (I4.16|) . we finally complete the proof. 


□ 


Lemma 4.5. 

Til (To) + 7l2(To) < CCo + CCo(C^^^ + Co4l2(To))jli(To)3/2 
+ CCoA,(Tof + CC^A,(Tof/^(alV^Bllo) 

+ C [ ° [ \Vu\^+ \\7B\^dxdt + C [ ° [ a{\\7uf+ \VB\‘^)dxdt. 
Jo Jr3 Jo Jr3 
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Proof. Combining Lemma 14.21 Lemma 14.31 and Lemma 14.41 we will have 


Ai(ro) + A 2 {To) < CCo + C f ° f a\Vuf dxdt + C f ° f |Vw|^ dxdt 

Jo JR3 Jo JR3 

(4.19) +C [ ° [ a{\u\^ + + \u\^ + |up(|VBp + \yu\^))dxdt 

Jo Jr^ 

+ C [ ” [ (|B|2 + |n|2)(|VBp + iVnp) dxdt. 

Jo Jr^ 


Next we need to estimate some typical terms on the right hand side. 

Terxnl: j;i'° dxdt 

n fTo / r \ 1/3 / /• \ 2/3 

\B\’^\y B\^ dxdt < / / \B\^dx] / \VB\^dx] dt 

;3 Jo \Jr3 j \Jr3 j 

n pTo !■ 

\VB\^dxdt + C / \Bf'dxdt 
;3 Jo JR3 

C f ° ( f \VBfdx') dt + C [ ° f \VBfdxdt 
Jo \dR3 / Jo dR3 

<c( f ivspdx) [" [ \VB\‘^dxdt + C [ ° f \VB\^dxdt 
\Jr^ j Jo Jr3 Jo Jr3 

< CCoAiiTof + C [ " [ \VB\^dxdt. 

Jo JR3 

Term 2 : J^° a\B\'^\Bt\'^ dxdt 


< 


< 


< 



rTo 
>0 V4r3 



\ 2/3 


\ 2/3 

cr|i?tp dxdt j 


<cCo/'^i(ro)4/3A2(ro)'/3. 


Term 3 : cr|i?p|Mp|Vilp dxdt 


a\B\‘^\u\‘^\yB\'^ dxdt < ^ ^ .^gdiS , i„,i8 


0 JR3 


< 


0 JR3 


a{\B\^+ \u\^+ \VB\*)dxdt 


r [ a\yB\Uxdt + {\\u\\l. + \\B\\U r<yi\\u\\U + \\B\\U)dt 
Jo JR3 JO 


Since 


\\u\\l. + ||i?||i4 < \\uU4^u\\l. + \\B\\^4yB\\l. 
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and 

rTo 

Jo 


a\\u\\U + ^\\B\\U dt < r a\\u\\U\Vu\\% + a\\B\\U\yB\\l, dt 

JO 

< [ \\\yu\\u\\puu^ + \\p{p)-pmLerdt+ [ \\\yB\\i4v^B\\i,dt 

Jo Jo 

pTo 


< a\\pi 


lL2 


IIVwlli^ dt + + CCoa\\V^B\\l,, 


we can get the estimate 

/■To /. 


r-To 


/ / a\B\'^\u\'^\VB\'Uxdt< I I a\VB\^dxdt 

Jo Jr3 Jq Jk3 

+CC'oAi(To)3/2(Co"/" + CoallV^SlIi, + CoA^iTo)). 

where we used Lemma 14.21 and the definition of Ai, A 2 . 

Using similar procedure as did for Term 1, we have 

/■To /■ 


(4.20) 


/o ./R3 


|U|" + |Mn(|VB|" + \Vundxdt 

fTo 


< CCoAiiTof + C [ " [ \VB\^dxdt + C [ " [ \\7u\^ dxdt. 
Jo Jr3 Jo 7r3 


Using similar procedure as did for Term 2, we have 

/■To 

(4.21) / 


a{\uf + \Bf){\Btf + \uf)dxdt < CCo/'Ai(ro)4/3A2(ro)i/^ 

Jo Jr3 

Using similar procedure as did for Term 3, we have 

/■To /■ 

|2 I I d|2\|„,|2/|v7 t}\2 


h Jkj 


(4.22) 


a{\uY + |U|")|u|"(|VU|" + |Vur) dxdt 
< CCoAi{T^f/\cl'^ + CoallVasili, + C'oA2(ro)) 
+ f I a\\7B\‘^ dxdt. 


Jo 4R3 

At last substitute (14.201) . (I4.21|) . (14.221) into (I4.19P and do simple reduction, we will 
complete our proof. □ 


Before going to the following part, we need to introduce some notations. Let 


(4.23) 

F = (A + ^)divM - P{p) + P{p), 

and 


(4.24) 

= <-<., j,fc = l,2,3. 


Through simple calculations, we could find that 

(4.25) = {pu%^ - {pu%. - {B ■ VB^ ),, + {B ■ WB %^, 
and 

(4.26) AF = div5 

where = pii^ + {\\B\’^)xj — div(B^B) with j = 1, 2, 3. 
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Lemma 4.6. 


E{Ta)<CCa + PE{A{To)), 

where ^(To) = Ai[Tq) + ^2(20) and Pe is a polynomial function with 2 as the 
lowest order. 

Proof. For the term a\\7B\‘^\u\‘^ dx, we have 

' |2 L. |2 ^ ^IL, II2 


(4.27) 


I a\VB\^\u\^ dx < a\\u\\i^ \VB\^dx 

R3 JR2 


<C'^i(To)(a||u||i4 + a||Vu||i4). 

For term cr||u ||^4 which appeared in the above inequality, we have 

(4.28) a||u||i4 < a\\u\\%^\\Vu\\%^ < A^{Tof/\ 

The analysis about cr||Vu ||^4 is a little bit more complex, by (I4.23|) . we know that 

(4.29) a||Vu||i4 < a||F||i4 + a||a;||i4 + a||P(p) - P(p)||i4. 

For the first term on the right hand side of the above inequality, we have 
cT||P||i4<a||F||^/"||VFf/." 

(4.30) < C(||Vu||l 2 + ||p-p||l2)1/" ( / a\u\'^ dx + f alWBl'^lBl^ dx 

\2r3 jh3 / 

< C (^Ai{To)^/^ + + F;(To) 3 / 4 ) . 

For the second term on the right hand side of (I4.29|l . we have 


3/4 


< crWuj 


,1/2 


0 ’|| W || 2,4 ^ 12 || U /||^2 

< llVuIlK" 


3/2 


l|Va;||^2 


(4.31) 


cr|Va;f dx I 


\ 3/4 


<Ai{ToY/^(I a\u\’^dx+ ! a\VB\^\B\^dx 
\Jw? 2r3 

< p1i(To)1/4 (A2(To)3/4 + F;(To)3/4) . 


3/4 


For the third term on the right hand side of (I4.29F we have 

(4.32) a\\P{p)-P{p)\\l.<c(^jjp-p\^dxy <CCl/\ 
Substitute (j4.30|l . (I4.31F (14.321) into (I4.29L we will obtain 

(4.33) cT||Vu||i4 < CCl'^ + C(kli(ro)'/" + Cy\A2{Tof/^ + E{Tof/^) 
Combining (|4.28F (14.3311 and ()4.27|1 . we finally have 


(4.34) 


a\^B\M^dx<CA^{n){cl/^A^{Tof/^+ C^ 

; \ 


1/2 

0 


(kli(To)'^^ + Cy\A 2 {nf/^ + E{nf/^)). 
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Using similar arguments as for term cr|Vi?p|Mp dx, we have 
I a\\/u\^\B\^ dx + I a\VBf\Bfdx 

(4.35) 

< CA,{To)(^Co^^A,iTof/^ + AiiToy/\A2{To)+EiTo)f/^). 
Summing up (14.341) and (14.351) . we obtain 
Ein) < Cyli(To)(Co^/"Ai(To)3/4 + 

+ (4ii(ro)i/4 + cd/")(4i2(ro)3/4 + i;(ro)3/4)) 

< Cl^^Aiiny/^ + Cl^^AiiTo) + Ai(To)5/4A2(To)3/4 
+ A,iTo)^/^E{Tof/^ + ^i(Ti)Cd/^Pl2(ro)3/4 + cl^Ui{To)EiTof/\ 
Using Young’s inequality, we can get 

, , E{To) <CC^Ai{To)i + CcjA,{To) + CCoAi{To)^ 

(4.36) j 

+ CAiiTof+CAi{To)^A2{To)i+CC~^AiiTo)A2iTo)"^. 

Through (I4.36L we could easily obtain our desired result. □ 

At this stage, we could give the estimate about term cr||V^i ?|||2 as follows 

'dx<C f a\Bt\^dx + C f cr\VB\^\u\^ + a\\7u\^\B\^ dx 
7r3 7r3 

<CA2(To)+CA(ro) 

< C(1 + Cd/')A2(To) + CAi(To)5/4yl2(To)3/4 

+ CCd/^Ai(To)A2(To)3/4 + CAi(ro)'^ + CCoAl(To)^ 
where we used Lemma 14.61 
Lemma 4.7. 

H{To) < CPnciCo) + CPha{A{To)), 

where Phc is a polynomial function with lowest order | and Pha is a polynomial 
function with lowest order |. 

Proof. We need to estimate every term appeared in the definition of H{To). For 
Io° /r 3 dxdt, we have 




/R3 


(4.37) 


a|VB|^ dxdt < C / a\\VB\\L2\\V^B\\l2 dt 


(4.38) 


< CA,{Toy/^a^/^V^BU2 ^ a\\V^B\\l2 dt. 

Jo 


For the last term of the above inequality, we have 

rTo 


(4.39) 


f ° f a\W^B\^dxdt<C f ° f aUBtl"^+ \WB\^\u\'^+ \Wu\'^\B\'^)dxdt 
Jo Jr3 Jo dR3 

<Ai(To)+ [ ° [ a{\yB\^\u\^ + \yu\^\B\^)dxdt. 
Jo Jr3 
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For the last two terms in the above inequality, we firstly have 

/•To /. 

cr|Vi?|^|u|^ dxdt 


/o JR3 
rTo 


< 


< 


\ 1/2 / r \ 

a\yB\'^dx] ij a\u\'^dx] dt 

\ / To ^ 

a\VB\'^dxdt\ if f alu]'^ dxdt 

'■ J \Jo j 


\ 1/2 


'R3 


fTo 


(4.40) 


1/2 


/ fTo / p \ 1/2 / /• \ 3/2 

<H{ToY/^\j \u\‘^dxj \Vu\^dxj dt 

( [ \yu\^dx 

\Jr3 


1/4 / /.To 


1/2 


1/2 


10 


|Vup dxdt 


Similarly, we could get 

/•To 

(4.41) 


a\\7u\^\B\^ dxdt < H{Toy/^Co^'^Ai{Toy/^. 

Jo Jr.3 

Substitute (j4.40|) and (14.4111 into (j4.39L we will have 
rTo 


(4.42) [° [ a\V^B\^dxdt<CAi{To) + CH{To)^^'^Co^‘^Ai{Toy/‘^. 

Jo Js.^ 

Substitute estimates (j4.37|) and (14.4211 into (I4.38|l . we obtain 

a\yB\^dxdt<CAi{To)^/^{{l + CoY/^A2iTo) 


10 Jr^ 


(4.43) 


+ A.iTof/^A^inf/^ + Co^^A,{To)A2iTof^^ 
+ A,iTof + CoA,iTo)y/^iA,iTo) 

+ Co'/%(ro)i/4i7(To)i/2) 


For the term crj Vu|"*’dxdt, we have 


(4.44) f [ (T|Vu|^dxdt< f f (t(|p — p|^ + |F|^ + |a;|"^) dxdt. 

Jo Jm.^ Jo Js3 

Concerning the first term on the right hand side of the above inequality, we have 

(4.45) f [ a\p - p\'^ dxdt < Co + f ( (j\F\'^ 

Jo JRJ Jo Jr? 


dxdt. 


10 Jr^ 

For the second term on the right hand side of (I4.44F we have 

rTo 


(4.46) 


/7 

10 dR3 


cr|J^|^dxdt< / °cr||F||i 2 ||VF||| 2 dt 
Jo 


< 


'R3 


\F\^dx] / cr|VFpdx / / iVF’pdxdt. 

/ \dR3 / Jo dR3 
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Due to 

(4.47) 

(4.48) 
and 

(4.49) 


'R3 


\F\^dx 


1/2 / „ \ 1/2 / „ 

<( / iVwpda;) +( / \p - p\^ dx 

\Jm? ) Vjr3 

<C(C'o + 41i(To))1/2, 

[ a\yF\^dx) ^ <C [ a{\u\^+ \yB\'^\B\^)dx 

JR3 / 7r3 

<CA2{To)+CE{To), 


\ 1/2 




f° f \VF\‘^dxdt<C f " f + \VB\‘^\B\‘^dxdt 
Jo Jr3 Jo JR3 

<CAi{Tq) + C ( f ° f \WB\^dxdt\ ([° [ \B\^dxdt^ 

yjo Jr3 J yjo JR3 j 


1/3 


< CAi(ro) + 

we can refine the estimate (14.4711 to 


r-To 


(4.50) 


/o ^R3 


a\F\Uxdt <C{Cq + Ai{To))^^HA2{To) + E{To)){Ai{To) 


+ i7(To)2/3Co/2Ai(To)i/®) 

Concerning the last term on the right hand side of (14.4411 . we have 

r'^o r /*Tb 

/ / CT|w|'*dxdt< / CT||a;||/,2||a;||^6 dt 

do dR3 do 

< ||a;||i 2 f cr|Va;pda; [ f iVwpdxdt 

dR3 do dR3 

< C24i(To)i/^(242(To)+£;(To))f [ “ [ p\u\^ dxdt 

^ Jo dR3 


(4.51) 


(•To 


+ 


/O ~'R3 


\VB\^\B\'^dxdt 


< 


CA,iTo)^/^{A2{To) + EiTo))(^Ai{To) 

+ (/'’/ |VB|3dxdt)^''^( [ ° [ \B\^dxdty^^) 

Jo dR3 do dR3 / 

< CA,{To)^/^{A2{To) + E{To))iAi{To) + Co^^H{Tof/^A,{Toy/^). 

Substitute (j4.45L (14.5011 and (14.5111 into (14.4411 . we obtain that 


rTo 


10 dR3 


(4.52) 


a\Wuf dxdt < CCo + C{Co + 24 i(To)) 1 / 2 (A 2 (To) 

+ EiTo))iAATo) + C^J^A,iToy/^H{Tof/^) 
+ CA,iToy/\A2iTo) + i?(ro))(24i(To) 

+ Co/'24i(To)1/3^(To)2/3). 
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Concerning the term f^° | V-Bp dxdt, we have 


rTo r rTi 

|3 


/o 


|VB|^ dxdt < C I II VBlI^i^ll dt 


rTo 


< 


(4.53) 


iVBlIi.dt) 

/ r'To 


1/4 


rTo 




V^BlIi.dt 


3/4 


^CIlVBlIi^ 


\ 1/4 

|VBp dxdt I 


10 




rTo r \ 

/ / IV^Bpdccdt) 

/o JK3 




< 


< 


CAi(ro)i/2Co/M f ° f \Bt\^ + \yB\^\u\^ + \yu\'^\B\^dxdt] 

\Jo Jr 3 j 

(Ai(To) 3/4 + Co/3Ai(ro)4/3iJ(To)i/3^ . 


\ 3/4 


For the term J^o |Vup dxdt, we have 


rTo 


|Vu|^ dxdt < 


dTo) 


10 


\Vu\^dxdt+ j f iVul^dccdt 

J(7{To) Jr^ 


r^{To) 


< 


(4.54) 


Vu\\%\\\pu\\%^ + ||P(p) - Pm%^) dt 


rTo 


la(To) JR^ 


\Vu\‘^dxdtj 


\ 1/2 


rTo 


la(To) JR^ 


{Vul'^dxd^ 


1/2 


< + 41i(To)3/2 + CCl^^H{To)^^^. 

Combining estimates (14.4311 . (I4.52F ()4.53ll . (14.5411 and (j4.36|l . we could finish the 
proof by a long but tedious calculations. □ 


Now, combining estimate (I4.37P , Lemma 14.61 Lemma 14.71 and Lemma 14.51 we 
will obtain 

(4.55) 7li(To) + A2iTo) < CPac{Co) + CPAA{Ai{n) + A^iTo)), 

where Pac{') is a polynomial function with lowest order | and Paa{') is a poly¬ 
nomial function with lowest order |. Since Ai{Trj) + A 2 {Tq) < and Co = 
WPo-RWh + WnoWlo + WBorHO < eo, if eo > 0 is small enough, we have 

(4.56) A,{To) + A2{To) < Cey\y^ + Cey^\y^ < 

It remains to prove the lower and upper bound of the density. Set F = log(p) 
which satisfies 


(A + p)T + [P{p) - P{p)) = -F. 
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For 0 < t < cr(To), we will have 


r<^(To) 


pg{Tq) 

\\F\\L^dt<C ||F||^/2||VF||y/di 

Jo 


(4.57) 


r<^iTo) 


< c / {\\v~pur^. +\\B- + \\b . m 

Jo 

pa‘{T'o ) 

Jo 


By estimate (14.561) . we have 


(4.58) 


r-^lTo) 


||Vpu||^/"||Vu||^/"d^<Cey^ 


r'^iTo) 


B-yB\\li^\\Vu\\]^2^dt 


r<yiTo) 


t-i/2(a||ij.Vi?||i.)'/"(a||Vw||i.)i/4rft 

(To) 


p(y{J-o) 

(4.59) <Cel^^ / f^/^((T||Vu||i2)^/^dt 

Jo 


<Ceo 


1/4 ( r^^"\-v^dt 

Jo t 


3/4 


r'yiTo) 


1/4 


cr||Vu|||2 dt 


<Cey^ey® < Ce^/®, 


(4.60) 


r-'^CTo) 


^(7(To) 
^0 




< 


r<^(To) 




r'^(T’o) 


®/® / r-{To) \ ®/® 




( r<^{To) \ / /-^(To) \ 

+ C'ei/ 4 (y t-^/^dt\ ij t\\VB\\l,dt\ 

^ 3/8 1/2 „ 1/4 9/16 ^ 3/4 

^ocg eQ +oeQ Cg ^ oeg > 
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and 


rcriTo) 




r<y(To) 


< 




(4.61) 


dTo) 


t-^/\t\\B-VB\\l.f/H^/^\\VB\\L^dt 


<cey^ 


/•o'(To) pcr^^o) 

/ t-^/^{t\\V‘^B\\l2f/^dt + Cel^^ t-^/V/^\\VB\\L4dt 

Jo Jo 


<^iTo) 


^ 3/8 3/16 ^1/4 9/16 ^9/16 

Soeg Cq +Oeo eg S oeg , 

where we used estimate (14.371) and l|4.56p . Substitute (|4.58p . (14.591) . (14.601) and 
(14.611) into (I4.57p . we obtain 

(4.62) / ”||F||io.d^<Cey^ 

Jo 

which implies that for t < a(To), 

(4.63) inf(logpo(a:)) - - Ct < logp{t,x) < sup(log(/Oo(a;))) + + Ct. 

So we can choose eo, t small enough such that for t < r < cr(Tg), 

(4.64) 


3 3 

-Co < p{t,x) < -CgV 


For T < ti < t 2 < Tq, we have 


/ ' ||F|U=o dt<c f ^Wy^puC^ + \\B ■ VB\\li^)i\\Vu\\li^ + \\B ■ VBWli^) dt 
Jti Jti 

< C{t2 - ti)i/2 r WVP^Wh + \\B ■ VSlIi. + llVulli, + \\B ■ VBlIie dt. 
Jti 

Now, we estimate the four integral terms on the right hand side of the above 
inequality. For the first two terms, we have 

(4.65) /'||Cpu||i.d^<7l2(Tg)<ey^ /hvu||i2dt<dl2(ro) 

Jti Jti 

The estimate about the third term appeared in (I4.49F so we have 

(4.66) r \\B-VB\\l2dt= r j \B\^\VB\'^dxdt<Cty‘^. 

Jti Jti 

For the fourth term, we have 


— ^0 ’ 


\B-VB\\l^ uu 


dt< f^\\VB\\%^\\V^B\\iydt+ f^WVBWydt 

Jti Jti 

V^BW^dt] 


\VB\\i2 


dt 


1/4 


t2 \ 3/4 


Ce, 


1/2 


<Cey\y^ + Cey^<Cey\ 


(4.67) 
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Combining estimates (I4.65L (|4.66|1 and (I4.67L we obtain 
f^jjFho^dt<c(t2-hy^Vo^^ 

Jti 


(4.68) 

< r]{t2 - h) + C^eo, 

where 77 > 0 is a small enough positive number. With estimate ()4.68p at hand, we 
can mimic the Navier-Stokes case [ 6 ] to obtain 


(4.69) 


3 3 

-Co < p{t,x) < -Co\ 


where 0 < t < Tq. With estimate (14.561) and (14.691) , we can complete the proof by 
continuity argument. □ 

5. Blow up Criterion and the Global Well-Posedness 

In this part, we give a blow up criterion and then prove the local solution can 
be extended to a global one. Firstly, let me give the blow up criterion as follows. 

Theorem 5.1. For dimension N = 3, let (p,u,B) be a solution of system 
satisfying 

p(0)>0, p-p€C{[0,T];H^), 

u,B G C{[0,T];H^)nL‘^{0,T-, H^). 

Let T* he the maximal existence time of the solution. If T* < +c», then it is 
necessary that 

limsup (||p(t)||L~ + \\u{t)\\Li + \\B{t)\\Li) = +00, 

tfT* 

for any q > 6. 

Proof. We use the contradiction argument. Assume that T* < +c» and 

(5.1) sup {\\p(t)\\L<=^ + \\u{t)\\Li + \\B{t)\\Li) = M <+ 00 . 
te[o,T*) 

In what follows, we denote C to be a constant depending on T, M, |luo||_y 2 , ||i?o||H 2 , 
\\po — p\\ff 2 . Firstly, from the energy estimates, we have 

(5.2) [ \p - p\'^ + p\u\'^ + \B\'^ dx + [ [ \yu\^ + \yB\^dxdt<C. 

JR3 Jo dR3 

Considering both eu and (Ol) . for r S [ 2 , c»], we have 


(5.3) 


Up - pIIl“(l'-) + IIv^“IIl“(l 2) + ||i?|| 


L“(L 2 


+ + I|v^IIl|(l2) < c. 

Let V = L~^VP{p) to be a solution of the following elliptic system 

(5.4) Lv := pAv + AVdivu = VP{p). 

By elliptic estimate, for r G [2, c»], we can obtain 

(5.5) \\VvU^<C\\P{p)-Pip)\\L^<C, \\V\\\L^<C\\yp\\L^. 

Now, we introduce a new unknown w = u — v. We can easily know that 

(A + /i)Adivi(; = (A + /i)Au - A{P{p) - P{p)), 
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and divw is called effective viscous flux which was used by many authorsigilllllSl 
[16] . Through simple calculations, we have 

(5.6) pdtw - pAw - AVdivw = pF + B - VB - iv(|Bp), 
where 

(5.7) F = -u-Vu + L-^Vdiv{P{p)u) + L"V((pP'(p) - P(p))divu). 
Multiply (15.6|) with dtw and integrating by parts, we have 

dt I /.tlVwp + AldivwP da: + - [ pldtwl'^ dx 

(5.8) Jr^ 2y„3 ' 

<C\\^F\\l,+C\\B-VB\\l,. 

For the first term on the right hand side of the above inequality, we can estimate 
as in the Navier-Stokes equations [SI 131] to obtain 

(5.9) \\^F\\l2 < (7(1 + ||V'u||i2) + e|l^atu;||i2, 

where e is a small enough positive number. Next, we need to estimate \\B ■ Vi?||i 2 , 
\\u-\7B\\i^2, ||i?-Vu||i 2 and \\u-Vu\\l 2 to make our later estimate more clear. For 
the term \\B ■ Vi?||i 2 , we have 

||B-VB||i2 <C||B||L,||VB||i2,/(,-2) 

<C\\yBU2+e\\y^BU2 

< C\\VB\\l 2 + e\\dtB\\L 2 + e||u • VBIU 2 + e\\B ■ Vu\\l 2 . 


Similarly, we have 

||u • VB\\l 2 < C|| VB||z ,2 + e\\dtB\\L 2 + e\\u • VB\\l 2 + e\\B ■ Vu\\l 2 . 

As in the Navier-Stokes equations [311131], we can easily obtain 
\\B ■ Vu\\l2 <C + C||Vu|U2 + e\\,/^dtw\\L2 
\\u ■ Vu||i 2 <C + C\\Vu\\l 2 + e||^5tw;||L2. 

So summing up the above four estimates, we obtain 

\\B-VB\\l2 + \\B-Vu\\l2 + \\u-VB\\l2 + \\u-Vu\\l2 

< C + (7||VM||i2 + (7||Vi3||2,2 + e||9ti?||i2 + e\\y/pdtw\\]:^2. 

Substitute (15.9p and (15.101) into (15.8p . we will have 

dt / ulVrcp + A|divu>p da; + / pldtwl'^ dx 

(5.11) Jr3 7r3 

< C(1 + IIVu||i 2 + IIVBlIi^) + eW^dtwWl. + e\\dtB\\l.. 

Multiply dtB to the third equation of system (HU and integrate by parts, we will 
obtain 

dt [ iVBl^dx+Z \dtB\^dx 

(5.12) J^3 

< C{1 + WVuWl. + WVBWl.) + eW^dtwWl, + e\\dtB\\l 2 , 
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where we used estimate ()5.10|1 . If e is small enough, summing up (j5.11|) and (j5.12|) . 
we will have 

dt [ + A|divuip + |Vi?p dx + f p\dtw\^ dx 

Jr3 

+ f \dtB\^dx<Cil + \\Vu\\h + \\VB\\l,). 

JR 3 

Integrate the above inequality for time variable, we will obtain 

+\\dtB\\LUL^) + \\^M\LUL^)<C. 

The above inequality combined with (15.51) implies that 

(5.14) + ||Vu||j;,|(j;,r) < C for re [2,6]. 

For the term || V^i3||2,2,(L2), we will have 

< C||ati3L|(i2) + ||u • Vi?L^(^2) + ||i? • VuL|(,,2) 


(5.15) 


<C< +00, 


where we used (|5.10|) and ()5.13p . Hence, we have ||VH||j;^|(^ 2 ) < C < +oo by 
interpolation. Now let us turn to the high order energy estimate. From the 
energy estimate iLemma l4.4F we have 


a(t){p\u\’^ + \Bt\^) dx + 


/R3 


<CCo + C 


cr(t)(|V'up + \WBt\^)dxdt 


(5.16) 


/O JR3 


cr(t)|Vu|‘* dxdt + C 


C 


/o JW3 

|2| d|2n 


|Vm|^ dxdt 


/o JR3 

i-T 


c 


pAw + AVdivic = pii — B ■ VB + -V(|i?p), 


(|VH|2|Hp + |VH|2 |u| 2 + |Vu|^|B|^) dxdt 
cr(t)(|up + |Hp)(|i?ip + [up + jupdVBp + jVup)) dxdt. 

Jo JR3 

Noting that 

(5.17) 

elliptic estimate yields that 

(5.18) \\^M\l3 < C\\pu\\l2 + C\\B ■ VB\\l2. 

From (15.101) and (|5.13|) . we obtain 

(5.19) f [ \yB\^\B\^ + \VB\^\u\^ + \yu\^\B\^dxdt<C<+oo. 

Jo JR3 
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Since 


/o Js.3 


< 


tT(t)|i?p|up|Vi?p dxdt 


(5.20) 


iBl^M^dx 


'R3 




WBfdx 


1/3 


dt 


< sup 

0<t<T 


\ 2/3 nT 

\B\^\ufdx\ I \\VB\\l<,dt 


< < C < +00 


we can use similar methods to obtain 

pT 


(5.21) [ [ a{t){\u\‘^+ \Bf)\u\‘^{\VB\‘^+ \Wuf)dxdt<C <+oo. 

Jo JR3 

Similar to the Navier-Stokes system [6], we can obtain 

(5.22) ||Vu||i4 < C||Vw||ie(l + IIVp^IUO, llVullia < c\\yu\\%\ 

Plugging ()5.19|) . (15.211) and (15.221) into ()5.16l) and noting ||Vu(t)||i6 G L^(0,T) by 
(I5.14P and (I5.15L we deduce by Gronwall’s inequality that 

(5.23) [ cr{t)ip\uf + \Btf) dx + [ [ ait){\Vuf+ \yBtf)dxdt<C. 

7r3 Jo jr3 

From the above inequality, elliptic estimate (15.1811 and Sobolev inequality, we have 

(5.24) l|V^u;||i 2 ^(ir) < C, for rG [2,6]. 

Using same methods as for Navier-Stokes system [B], we have 


(5.25) 

Since 


and 


< C, \\Vp\\Lr<C, for rG[2,6]. 


< CiWBth^ + \\u\\l^\\VB\\l^ + ||i?||Loo|lVu|U2) 

<c + l\\yML^ + \\\v^BU2, 


\^^u\\L2<c\\y^w\\L2 + c\\y\\\L- 

<C(||ph|U2 + ||il|Uo=|lVi?||i2 + ||Vp|U2) 

<c + 1\\v^b\\l2, 


we have 
(5.26) 


\\V^u\\l2 + \\V^BU2<C. 

By the above estimates, we easily know ||w||lo° -|- ||i?||L°° < G by interpolation, so 
we will obtain 

+ ||B|U^||V2u|L^(i2) + ||VnVi3L|(i2)) < C, 


(5.27) 
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where we also used 


X 1/2 \ 1/2 


< C||V«|U^(i4)||V«L|(i4) + C||VS|U^(^4)||VSL^(i4) 
<C< +c». 


Then, similar to Navier-Stokes system, we can get 

/•T 

(5.28) 


f W^^Mlh dt < C. 
Jo 


From (I5.5F (I5.26L ()5.27p and (I5.28L we will obtain 

l|vV/)lli= + I|v^i3(t)lli. + r \\vMh + Iiv^i3||i. dt 

(5.29) ^ 

<C + C f \\V^p{t)\\l,dt. 

Jo 

Form the continuity equation, we have 

(5 30) + ||Vu(t)||L~)||V2(p(t) - p)\\h 

+ c\\vMmh. 

Summing up (15.2911 and (Ih.dOL we conclude by Gronwall’s inequality that for 0 < 
t < T\ 

l|p(/) - + l|ii(/)lli?2 + \\B{t)\\H'^ < C. 


This ensures that the solution can be continued after t = T*. □ 

After we get the blow up criterion, we can give the proof about Theorem 12.81 
From Theorem 12.61 Proposition 13.31 and Proposition 13.41 we obtain a solution 
{p, u, B) of (11.11) satisfying 

y < P < 2co-i, p-pe L’^{Bl[? n ii^), 

u,bg n n 4(5^/^' n ). 

Moreover, it holds that 

\\P ~ P\\— *//(IIPo — + l|wo|lij-« + II-5 o||^-5)i 

||(m, i?)||£^(g-S) + \\{u,B)\\i^^g2-s-^ 

< C{1 + ||ao||^3/p)(||Mo|lij-i + \\Bo\\h-s + TP+||po - pllffi-?), 

P.l 


(5.31) 

(5.32) 


(5.33) 

(5.34) 

(5.35) 


^ 2||ao||^3/p + 

^P,l 
C2 + C3 


I"^11 r oof r3/p- 1\ H" ll'^li f 1 ( D3/P+l^ ^ 




- (1 + 


2c3 


Ikoll 


d3/p J 

•^P.l 
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Here and in what follows, the constant C depends only on A,/i, co, p,p, s. Recall 
that the initial data satisfies 

C2 II n II ^ C2 


Ipo - pIIl2 < Cl, ||UollzS3/p< 


B, 


C.”" (l + llaolUa/.)^’ (l + ||ao|U3/.)2’ 


||'*4o||_h--« < 


C2 


(1 + ||ao||^3/.)^(l + IIpo - prH2 + \\uo\\%^ + ||i?o|II4") 


,4/3. 


= C2, 




C3 


(1 + llaoll33/.)V3(1 + IIpo - + \\uo\\%^ + \\Bo\\%^ 


= C3- 


Taking C 2 small enough, by Remark ?? and (I3.24I1 . the existence time T has a lower 
bound 

c 


(5.36) 


T > 


(l + ||ao||^3/.)4(l + ||ao||^3)i2’ 

where c is a small positive number. 

Now, for any Ti < T, we have 

II '®)IIl“ (B'j) + II 

< 1^(1 + ||ao||g3/p)(c2 + C3 + P+Ti||po — Pllffi—5); 
by the conditions on the solution just mentioned. For r G (1,2 — ^), we have 

+ CXD 


(5.37) 


(5.38) 


< ^2^Ti/n|A,u|U~(L3) +E2^II^'«II™)II^'«I^ 




<c^(hiir,^(S2-:2^) + ihiih^^(4V)) 

< C{1 + ||ao|| A 3 /p)(c 4 + C 5 + P+Ti||po — PII// 2 ). 

■°P.l 

Similarly, we have 

(5.39) I|w||l5,^(l2) < C{1 + llaoll b3/p)(c 4 + cg + T+Ti||po - pWh^), 

(5.40) l|VR||Lr^(L2) < (7(1 + ||ao||^3/p)(c4 + cg + P+TiUpo - p\\m)^ 
and 

(5.41) II^IIl5.^(l2) < (7(1 + llaoll^3/p)(c4 + cg + P+Ti||po - p\\m)- 
Then taking c = |, we have 

(5.42) ll(w> ^)ll4/"(ffi) ^ C'(l + ll«o||4/p)(c4 + eg + P+TiIIpo - p| 42 ). 
Hence, there exist to G {0,Ti) such that 


(5.43) ||(a,i?)(to)||4/i <C(l + ||ao||^3/p) + P+4'/^||po - p||ff 2 j . 

For density, using simply energy estimates, we will get 

(5.44) ||p(t) - p||i ,2 < C(1 + P+)(ci + 4'/"(||no||L2 + ||Po||l 2 )). 
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Take Ti as 
Ti = 


(5C)3(1 + P+)3(l + IIpo - + IIK, Bo)|li.)(l + llaoll’ 


such that Ti < T and 


C{l + \\a4^v,)P+Tl/^\\p^--p\\H^ < 

Pil 

c{i + p+)tI/\\\u4l^ + \\b4l^)<- 


eo 


Then, we choose ci, C 2 and C 3 small enough so that 

1^(1 + llooll r3/p)c2 1^(1 + ||ao|l 33 /p)C 3 


C(l + P+)ci < 

5 


T, 


2/3 


5 


P 


2/3 


-1 1 

Through this choice of ci, C 2 and C 3 , by (|5.43|1 and (I5.44L we have 

l|a(^o)||L2 + ||u(to)||ffi + \\B{to)\\H^ < eo- 
Theorem 14.11 implies that 


- 5 


Co 


< p < 2Co \ ||u(t)||L-(L6) + ||u(t)||L~(L6) < C. 


So the solution can be extended to a global one by Theorem 15.II 

6. Appendix 

In this section, for the reader’s convenience, we firstly give the Lagrangian trans¬ 
form for terms appearing in the MHD system as follows 

= j-idivj,(adj(DA)|Bn, 

B ■ VxB = diVx{B 0 P) = J“3clivy(adj(DA)^P 0 B), 

(6.1) div3,uP = J“^divy(adj(DA)ft)P, 

(divxM)H -I- u ■ VxB = Vx{uB) = J“^divy(adj(DX){tP), 

B ■ VxU — diYx{B 0 m) = J“^divy(adj(DA)^'u 0 B). 

Then, let us look at the following Lame system with nonconstant coefficients 

(6.2) dtv — 2adw{pD{v)) — 6V(Adivu) = /. 

We assume that the following uniform ellipticity condition is satisfied: 

(6.3) a := min 


inf (au)(t,x), inf (2an-|- 6 A)(t, a:) >0 

(i.3^)e[0.T]xR" (t.x)G[0,T]xR«‘ J 

Concerning this equation, we have the following lemma. 

Lemma 6.1. [IHIISS] Let a, b, A and p he bounded functions satisfies 16'. 51) . Assume 
that aVp, bVX, pVa and XVb are in L°°{0,T; Bp^^) for some 1 < p < 00 . There 
exist two constants p and k such that if for some m € Z we have 


inf Smifiap + bX){t,x), 


inf Sm{ap){t,x) ) > 


(t,a:)G[0,T]xR'^ ■ .(t,x)G[0,T] xR'^ 

||(/ - Sm){pVa,aVp,XVb,bVX)\\^^^^^/-p-i^ < pa 


2 ’ 


(6.4) 
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then the solutions to lih.Si) satisfy for all t € [0,T], 

< C (||z;o|Ib| + II/IIl1(b|)) exp ^ ||S'^(/4Va, aV^i, AV6,&VA)|||„/p , 

whenever — min(n/p, n/p') < s < n/p — 1 is satisfied. 

Lemma 6.2. [33] Let p G [1, +oo], 1 < P 2 < Pi < +oo, and let v solve 


I dtv - vAv = /, 

= vq- 

Denote p'^ = (1 + 1/pi — l/p 2 )“^. Then there exist two positive constants c and C 
depending only on N and such that 




'l _ g-c!vT2''''pi A 


' (jez 


ci/pi 


^2«(^-2+2/'’^)||AJ||,p. 


1 — e 


-cyT2^''pi \ 


qJ\\D^^(Lp) 


cvp[ 


In particular, we have 


- :yi/pi ' j^i/p'i iilJ''(B; 


C 


pa + 


c 


P2 ( oS-2+2/P2 \ ■ 


Moreover u belongs to C{[0,T]; B^). 
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